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ABSTRACT 

HD 156846b is a Jovian planet in a highly eccentric orbit (e = 0.85) with a period of 359.55 days. The 
pericenter passage at a distance of 0.16 AU is nearly aligned to our line of sight, offering an enhanced 
transit probability of 5.4% and a potentially rich probe of the dynamics of a cool planetary atmosphere 
impulsively heated during close approach to a bright star (V — 6.5). We present new radial velocity 
(RV) and photometric measurements of this star as part of the Transit Ephemeris Refinement and 
Monitoring Survey (TERMS). The RV measurements from Keck- HIRES reduce the predicted transit 
time uncertainty to 20 minutes, an order of magnitude improvement over the ephemeris from the 
discovery paper. We photometrically monitored a predicted transit window under relatively poor 
photometric conditions, from which our non-detection does not rule out a transiting geometry. We 
also present photometry that demonstrates stability at the millimag level over its rotational timescale. 
Subject headings: planetary systems - techniques: photometric - techniques: radial velocities - stars: 
individual (HD 156846) 



1. INTRODUCTION 

The discovery of exoplanets using the transit technique 
is becoming increasingly dominant amongst the various 
detection methods. Examples of major contributors to 
the ground-based discovery of transiting exoplanets are 
the Hungaria n Auto mated Telescope Net work (HATNet) 
(iBakos et alj I2004D and SuperWASP (|Pollacco et all 
2006). From the vantage-point of s pace, the major con- 
tributors are the Keple r mission (IBorucki et all 120101) 
and the CoRoT mission (jBarge et al.ll2008D . Thediscov- 
eries provided by these surveys are producing insights 
into the exoplanet mass-radius rel ationship, extendin g 
down towards super-Earth planets fScagcr ct al. 2007). 
Although these space-based surveys are expected to ex- 
tend the period sen sitivity to longer p eriods, such as the 
case of CoRoT-9b (jDeeg et all llOlO). the picture is in- 
complete since the surveys are strongly biased towards 
short-period planets around relatively faint host stars. 
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Several planets discovered with the radial velocity tech- 
nique have subsequently been found to transit, the first 
of which was HP 209458b (jCharbonneau et all l2000t 
iHenrv et al.ll2006T ) . The brightness of their host stars has 
facilitated further characterization of their atmospheres, 
such as the cases of HD 189733b and HD 149026b 
(e.g., iKnutson et all (l2009allb|) see also review arti- 
cle by iSeager fe Demind (J2010)) . The Neptune-mass 
planet orbiting GJ 436 became the fi rst known transitin g 
planet around an M-dwarf primary dGillon et al.1 120071 
The detection of tra nsits for the plan ets HD 17156b 
(iBarbieri et all 12001 and HD 80606b (jLaughlin et all 
1200ft iMoutou et aLll2009l). enabled by th eir high eccen- 
tricities (jKane fc von Braunl 120081 |2QQ9l), provided the 
first insights into the structures of longer-period plan- 
ets. Many of the known planets with orbital peri- 
ods larger than a few days have yet to be photomet- 
rically monitored at predicted transit times, hampered 
mostly by insufficient orbital parameter precision to ac- 
curately predict when the planet might transit. Fur- 
ther discoveries of long-period planetary transits around 
bright stars are vital to understanding the dependence 
of planetary structure and atmosphe ric dynamics on the 
periastron distance of the planet dFortney et all 120101 : 
iKane fc Gelinol l2010t iLangton fc Laughlinl \20m . Pro- 
vided the orbital parameters can be determined with suf- 
ficient precision, monitoring planets detected via the ra- 
dial velocity technique at predicted transit times provides 
a means to increase the sample of long-period transit- 
ing planets (|Kane et alll200ft 120101 ) . There exist efforts 
to detect transits of the known radial velocity planets, 
such as the Spitzer search for transits of low-mass planets 
(jGillon et al.ll2010f ). The Transit Ephemeris Refinement 
and Monitoring Survey (TERMS) is a program which is 
capable of monitoring long-period as well as short-period 
planets by refining the orbital parameters of the system. 

Here we present a detailed analysis of one such system. 
The massi ve planet orbit i ng th e star HD 156846 was dis- 
covered bv iTamuz et all (|2008f ) using the CORA1IE in- 
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strument. The planet is in a highly-eccentric orbit with 
a period of slightly less than a year. The periastron ar- 
gument of the orbit is such that the transit probability 
is significantly enhanced compared to an equivalent cir- 
cular orbit (5.4% compared to 0.9%). Our combined fit 
to new Keck data along with the discovery CORALIE 
data greatly improves the orbital parameters for the sys- 
tem, allowing an accurate prediction of possible transit 
times. We also find no evidence for additional compan- 
ions in the system through high-precision radial velocity 
(RV) data acquired during periastron passage over mul- 
tiple orbits. The long-term photometry presented here 
establishes the photometric stability of the host star. We 
present photometry acquired during a predicted transit 
window which places an upper-limit on a transit for this 
planet. Finally, we discuss additional constraints on the 
mass and orbit of the planet from a potential transit 
null-result. 

2. SCIENCE MOTIVATION 

Here we describe why the planet orbiting HD 156846 is 
a particularly interesting target and the potential gains 
which may be achieved through further studies. 

HD 156846 is an extraordinarily bright star (V = 6.5), 
brighter indeed by a factor of ~ 2.9 than either of the 
planet hosting stars HD 209458 and HD 189733. The op- 
portunities for follow-up studies of a fundamentally new 
type of planetary atmosphere would therefore be close 
to optimal. Note that massive, relatively cold planets 
such as this one have intrinsically difficult atmospheres to 
study via transmission spectroscopy. Their atmospheric 
scale heights are of order a factor of 20 smaller than 
typica l hot Jupiters (see for example IVidal-Madiar et aT] 
(|201 If ) ) . such that a bright host star is needed to achieve 
adequate signal-to-noise. 

Given the properties of this star (see Section l3~2|) . the 
received flux of the planet at apastron will be nearly iden- 
tical to the flux received by the Earth from the Sun. It 
is therefore not unreasonable to expect that the planet 
during this phase of the orbit will be sheathed in re- 
flective white water clouds. At some point prior to pe- 
riastron, when the received flux increases briefly to a 
value nearly 150 times that at apastron, the received flux 
should be sufficient t o flash the water clouds to steam 
(|Sudarskv et al.ll2005l ). Planets in the post- water cloud 
temperature regime are expected to have atmospheres 
transparent down to large pressure depths, which will 
cause a dramatic drop in the planet- wide albedo. Be- 
cause of the extra-bright primary, HD 156846 will al- 
ways represent one of the very best targets in the sky for 
reflected light observations of this time-sensitive albedo 
change, which gives real insight into the atmospheric dy- 
namics. This system is thus likely to gain in importance 
as photometric sensitivities improve, and so any knowl- 
edge of the inclination (and whether it transits) is very 
i mportant. 

Hamilto n fc Burns! (|1992f ) have shown that the Hill ra- 
dius at periastron is a good representation of the sta- 
bility zone for a satellite of a planet in very eccen- 
tric orbit around a star. The relatively large mass of 
HD 156846b leads to a sizable effective Hill Sphere at 
periastron (~ 0.02 AU, ~ 47 Jupiter radii), indicating 
that detectable Earth-mass satellites orbiting the planet 
can be dynamically stable over the lifetime of the system 



if the tidal quality factor, Q, is of o rder the Jovian value 
or higher (iBarnes fe O'B rien 2002). Given the mass ra- 
tios observed for the Jovian planets and their satellites 
in our own solar system, one might reasonably expect 
a ~ 0.5Af ffi satellite, which wo uld be readily detectable 
using transit timing techniques ( Kipping 2009) , and per- 
haps even directly via space-based photometry. Given 
the inflated transit probability for this planet, and its 
potentially interesting dynamical history, this becomes a 
prime candidate in this regard. 

3. KECK MEASUREMENTS AND REVISED ORBITAL 
PARAMETERS 

3.1. Observations 

We obs erved HD 156846 with the HIRES echelle spec- 
trometer ([Vogt et al.lll994f ) on the 10-m Keck I telescope 
with the goal of improving the accuracy of the predicted 
transit time to guide and temporarily anchor a photomet- 
ric monitoring campaign. Our Kec k observations post- 
date the CORALIE measurements (jTamuz et al.l 120081) 
and span 2009 May to 2010 October. The 41 Keck 
RV measurements were made from observations with an 
iodine cell mounted directly in front of the spectrom- 
eter entrance slit. The dense set of molecular absorp- 
tion lines imprinted on the stellar spectra provide a ro- 
bust wavelength fiducial against which Doppler shifts arc 
measured, as well as strong constraints on the shape of 
the spectro meter instrumental profile at the time of eac h 
observation (]Marcv fe Butlerlll992t [Valenti et alJll995f) . 
We measured the Doppler shift of each star-times-iodine 
spectrum us i ng a modelling proc edure descended from 
iButler et all (jl996l) as described in lHoward et~aT1 (|2009ft . 
The times of observation (in barycentric Julian days), 
relative RVs, and associated errors (excluding jitter) are 
listed in Table [T] In cases when we observed the star 
3-5 times in quick succession, we report the mean RV 
and appropriately reduced uncertainty. We also observed 
HD 156846 with the iodine cell removed to construct a 
stellar template spectrum for Doppler modelling and to 
derive stellar properties. 

3.2. Stellar Properties 

We used Spectroscopy Made Easy (IValenti fc Piskunovl 
1996) to fit high-resolution Keck-HIRES spectra of 
HD 156846 (HIP 84856, TYC 6242-00339-1), applying 
th e wavelength inter v als, l ine data, and methodology 
of IValenti fc Fischer! (|2005). We further constrained 
the surface g ravity using Yonsei-Ya le (Y 2 ) stellar struc- 
ture models (iDemarque et al.|l2004fl and revised Hippar- 
cos paralla xes (Ivan LeeuwenT I2007D . with the iterative 
method of IValenti et al.1 ([2009). The resulting stellar 
parameters listed in Table [2] are effective temperature, 
surface gravity, iron abundance, projected rotational ve- 
locity, mass, and radius. HD 156846 lies 1.13 mag above 
the Hipparc os average mai n sequence (My versus B — V) 
as defined bv lWrightl(|2005t ) . These properties are consis- 
tent with a metal-rich GO star evolved slightly off of the 
main sequence. The stellar radius, = 2.12±0.12 Rq, 
is crucial for estimating the depth and duration of a plan- 
etary transit. 

Our characterization is mostly consis t ent w ith the stel- 
lar properties reported bv lTamuz et~all([2008l ). Although 
small, key differences are the larger mass (M+ = 1.43M ) 
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TABLE 1 
Keck Radial Velocities 



Radial Velocity Uncertainty 



BJD - 2440000 


(ms- 1 ) 


(ms" 1 ) 


14957.02010 


-111.41 


1.94 


14983.95333 


-79.61 


1.73 


14984.87930 


-76.60 


1.74 


14985.97062 


-77.24 


1.63 


14986.96751 


-78.81 


1.69 


14987.91548 


-80.71 


1.54 


15015.85870 


-19.24 


1.60 


15016.78653 


-11.02 


1.61 


15026.95355 


22.14 


1.29 


15028.96769 


24.15 


1.27 


15041.87653 


81.78 


1.64 


15042.95724 


91.65 


1.63 


15043.79662 


94.00 


1.56 


15048.77399 


124.15 


1.90 


15073.72838 


561.70 


1.38 


15074.72238 


562.49 


1.22 


15075.72238 


507.32 


0.72 


15076.73007 


388.21 


0.78 


15077.72296 


222.03 


0.73 


15078.72435 


67.34 


0.59 


15078.75583 


61.16 


0.57 


15079.72234 


-60.59 


0.74 


15080.72602 


-144.09 


0.49 


15081.72293 


-211.59 


1.22 


15082.71457 


-252.11 


0.86 


15083.71752 


-278.63 


0.67 


15084.71590 


-308.49 


1.13 


15111.70803 


-337.84 


1.80 


15229.17709 


-213.95 


1.67 


1 KQon none? 
1532U. 98063 


-107.50 


1 a'? 
1.0 1 


15375.84930 


-20.06 


1.82 


15402.81296 


92.47 


1.65 


15410.78620 


144.58 


1.51 


15433.71936 


569.13 


1.20 


15435.71871 


457.38 


0.69 


15439.71931 


-97.88 


0.75 


15440.75164 


-174.32 


0.70 


15455.70798 


-363.19 


0.72 


15468.74647 


-344.65 


1.29 


15471.70625 


-340.87 


0.77 


15486.68538 


-316.77 


0.92 


TABLE 2 



Stellar. Properties 



Parameter 



Value 



M v 

B-V 

V 

Distance (pc) 
T cff (K) 

log 3 
[Fe/H] 

vsrai (kms x ) 
M* (M ) 
R* (R@) 
logR' aK 
Skk 



3.055 
0.557 
6.506 
49.0 ±2.2 
5969 ± 44 
3.92 ± 0.08 
+0.17 ±0.04 
5.05 ± 0.50 
1.35 ± 0.045 
2.12 ± 0.12 
-5.071 
0.144 



and higher ef f ective temperature (T c s — 6138 ± 36 K) in 
iTamuz et all (|2008D . In addition, we measured the stel- 
lar activity by measuring the strength of the Ca II H & 
K lines which give calibrated 5*hk values on the Mt. Wi l- 
son scale and log i?^ K values (jlsaacson fc Fischenl201(fl . 
The median of log i?y K and Shk values are listed in Table 
[2] and demonstrate the HD 156846 is chromospherically 



> 




2004 2006 2008 

Time (Years) 



2010 



E 



o 

CD 
> 



600 



400 



200 - 



-200 




4 6 
BJD - 2455073 




0.6 
2455078 



Fig. 1. — Time series RV measurements of HD 156846 from 
CORALIE (red filled circles) and Keck-HIRES (blue filled circles). 
The adopted Keplerian orbital solution (Table [3j is shown as a 
dashed line in all panels. Dotted vertical lines illustrate the times 
of ingress, mid-transit, and egress for a predicted time of transit 
on 2009 September 3 (see text). Top: all CORALIE and Keck- 
HIRES measurements. Middle: Keck measurements during the 
2009 September periastron passage. Bottom: Keck measurements 
on 2009 September 3. 

quiet, which is consistent with the photometric stability 
described in Section [5] 

3.3. Keplerian Models 

With the goal of improving the orbital solution for 
HD 156846b, we constructed single-planet Keplerian 
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models of the RVs using t h e orb it fitting techniques 
described in iHoward et al.l (|20100 and the partially 
linearized, least-squares fitting procedure described in 
iWright fc Howard (2009). Each velocity measurement 
was assigned a weight, vo, constructed from the quadra- 
ture sum of the measurement uncertainty (orv) an d a 



HD 156846 b 



P = 359.55 days 
a = 1.10 AU 
e = 0.85 
a = 51.6° 



jitter term (emitter), i-e. 



l/(a 



2 

RV 



a 



jitter/ 



We chose 



jitter values of critter = 3.38 and 5.90 ms 1 for Keck 
and CORALIE to satisfy the condition xt — 1 f° r Ke- 
plerian fits to those data sets individually. These val- 
ues are consistent with the expected jitter of a slightly 
evolved early G star observed with those instruments. 
Sources of jitter include stellar pulsation, magnetic ac- 
tivity, granulation, undetected planets, and ins trumental 
effects ((Isaacson fc Fische3l20Tol:IWrightll2005( ). 

The Keplerian parameter uncertainties fo r each planet 
were derived using a Monte Carlo method (|Marcv et al.l 
120051 ) and account for correlations between parameter 
errors. Specifically, our method accounts for correlations 
between T c and the other Keplerian parameters (notably 
e) to provide an accurate estimate of the transit time. 
Uncertainties in M sin i and a reflect uncertainties in M* 
and the orbital parameters. 
We considered models base d on three data sets: 

, [2010h and Keck-HIRES 
Each model con- 



CORALIE (jTamuz et al.l I200S 
(Table Q]) individually and combined, 
sists of a single planet in Keplerian motion with the 
parameters listed in Table [3J We allowed for an arbi- 
trary RV offset in the CORALIE measurements at JD 
2,454,279, the time of an instrument upgrade. We also 
allowed for an RV offset between the Keck-HIRES and 
CORALIE measurements. Our models include a linear 
velocity trend. Because the inclusion of a trend does not 
lower xti the data do not provide evidence for a distant 
third body in the system. Nevertheless, our models in- 
clude the trend to provide sufficient model flexibility to 
achieve properly estimated parameter uncertainties. 

Table [3] lists the parameters of the three, nearly iden- 
tical o rbital models. Our refitting of the iTamuz et al.l 
( 2008) data have a slightly lower rms than they reported 
because we excluded measurements with uncertainties 
greater then three times the median. The higher preci- 
sion Keck measurements yield a substantial improvement 
in the estimated parameters. The uncertainty in the pre- 
dicted time of transit on 2009 September 3 is improved 
by an order of magnitude, from 0.227 days to 0.023 days. 
(This improvement is also due to the timing of the Keck 
measurements to coincide with that epoch.) We adopt 
the CORALIE+Keck model which has a period uncer- 
tainty of one part in 50,000 and a transit mid-point un- 
certainty of only 20 minutes. 

Figure [T] shows the CORALIE and Keck measure- 
ments in time series with our adopted model overlain. 
The middle panel shows Keck measurements during the 
2009 September periastron passage. These measure- 
ments and similar ones during the 2010 periastron pas- 
sage provide substantial leverage to constrain the orbital 
period and transit time. The bottom panel shows the 
Keck measurements on 2009 September 3, when a tran- 
sit was predicted. The dotted vertical lines show the 
predicted mid-transit time and the earliest ingress and 
latest egress times that are consistent with the adopted 
CORALIE+Keck orbital solution. These times were 




Y to observed 



Fig. 2.— The orbit of the planet orbiting HD 156846 (solid line) 
and the orbits of the Solar System planets for comparison (dashed 
lines) . 

computed from t c ± (2<7t e + idur/2), where t c and <7t c 
are the predicted transit time and uncertainty, and idur 
= 490 minutes is the duration of an equatorial transit 
of the star with R* = 2.12 ± 0.12 R Q . Our Keck mea- 
surements that night were taken as early as possible but 
still missed the transit window and do not provide a con- 
straint on_the_inchn^ionfrom the Rossiter-McLaughlin 
effect (IGaudi fc Windl20"ol . 

4. TRANSIT EPHEMERIS REFINEMENT 

As described by iKane &; von Braun ( 2008), the tran- 
sit probability of a planet is a strong function of both 
the eccentricity and the argument of periastron. In par- 
ticular, the transit probability is the strongest when the 
periastron passage occurs close to the star-observer line 
of sight, or where u = 90°. The orbit of HD 156846b, 
shown in Figure [H is well-suited for photometric follow- 
up at predicted transit time since the alignment of the 
major axis with the observer clearly leads to an enhanced 
transit probability. 

From the derived stellar properties shown in Table [5] 
and the planetary properties from Table [31 we derive 
a planetary radius of R v = 1-l^j us ing the methods 
described in iBodcnhci mer et all (|2003l ) . This produces 
a transit probability of 5.4% and a predicted transit 
depth of 3 millimags. The uncertainty in the stellar 
mass/radius and subsequent uncertainty in the planetary 
mass/radius have a minor effect on the estimated transit 
duration but in no way affects the predicted transit mid- 
points since t hese are derived from the orbital parameters 
(jKane et al] 12009). Based upon the revised orbital pa- 
rameters, the transit ephemeris has been calculated for 
the period 2009-2016 and is shown in Table [U 

As described in the previous section, the size of the 
transit window for 2009 was 0.37 days, or ~ 9 hours. 
The uncertainty on the transit mid-point is small, ~ 20 
minutes, so the transit window is largely dominated by 
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TABLE 3 

Keplerian Orbital Models 



Parameter CORALIE Keck CORALIE+Keck (adopted) 

P (days) 359.5400 ±0.1072 359.5517 ± 0.0270 359.5546 ± 0.0071 

T c a (JD - 2,440,000) 15078.459 ± 0.219 15078.470 ± 0.017 15078.483 ± 0.015 

T p b (JD - 2,440,000) 15076.699 ± 0.227 15076.671 ± 0.023 15076.686 ± 0.021 

e 0.84696 ± 0.00109 0.84828 ± 0.00061 0.84785 ± 0.00050 

K (ms- 1 ) 464.03 ±2.14 463.64 ±1.20 464.14 ±0.96 

u (deg) 52.51 ±0.39 51.42 ±0.18 51.62 ±0.16 

dv/dt (ms-'yr -1 ) 3.51 ± 1.41 0.86 ± 1.36 1.55 ± 0.88 

Msini (M Jup ) 10.60±0.29 10.54 ± 0.29 10.57±0.29 

a (AU) ' 1.096 ±0.021 1.096 ±0.021 1.096 ±0.021 

rms (ms" 1 ) 6.85 3.45 6.06 



a Time of transit. 

b Time of periastron passage. 



TABLE 4 

Refined transit ephemeris for HD 156846b. 



Beginning 


Mid-point 


End 


JD 


Date 


JD 


Date 


JD 


Date 


2455078.30 
2455437.85 
2455797.39 
2456156.94 
2456516.49 
2456876.04 
2457235.58 
2457595.13 


2009 09 03 19 09 

2010 08 29 08 17 

2011 08 23 21 25 

2012 08 17 10 34 

2013 08 11 23 42 

2014 08 06 12 50 

2015 08 01 01 59 

2016 07 25 15 07 


2455078.48 
2455438.04 
2455797.59 
2456157.15 
2456516.70 
2456876.26 
2457235.81 
2457595.37 


2009 09 03 23 35 

2010 08 29 12 54 

2011 08 24 02 12 

2012 08 17 15 31 

2013 08 12 04 49 

2014 08 06 18 08 

2015 08 01 07 27 

2016 07 25 20 45 


2455078.67 
2455438.23 
2455797.79 
2456157.35 
2456516.91 
2456876.48 
2457236.04 
2457595.60 


2009 09 04 04 01 

2010 08 29 17 30 

2011 08 24 06 59 

2012 08 17 20 28 

2013 08 12 09 57 

2014 08 06 23 26 

2015 08 01 12 55 

2016 07 26 02 23 



Note. — The columns indicate the beginning, mid-point, and end of the transit window in both Julian and calendar date. The calendar 
date is expressed in UT and includes the year, month, day, hour, and minute. 

the transit duration. The small uncertainty on the pe- 
riod ensures that the size of the transit window does not 
grow substantially with time. In 2016, for example, the 
transit window is only 2.5 hours longer than it was in 
2009. The transit duration makes it very difficult to at- 
tempt complete coverage of the transit window without 
a multi-longitudinal campaign and exceptional weather 
conditions. However, ground-based observations of ei- 
ther ingress or egress will be feasible from a given ob- 
serving location, provided photometric precision require- 
ments are met. 

5. PHOTOMETRIC STABILITY 

Here we describe photometry acquired outside of the 
transit window for the purposes of studying the stabil- 
ity of the star. The only source of time-series photom- 
etry for HP 156846 in the li terature is the Hipparcos 
catalog (perrvman et al.lll997l ). HD 156846 (HIP 84856) 
was measured 82 times over the three-year duration of 
the Hipparcos mission; its variability classification in the 
catalog is blank, indicating that the star "could not be 
classified as variable or constant" . The scatter (standard 
deviation) of the 82 observations was 0.007 mag. 

We acquired new photometry of HD 156846 with the 
T8 0.80 m automatic photometric telesco pe (APT) a t 
Fairborn Observatory in southern Arizona (|Henrvi ll999). 
T8 uses a two-channel precision photometer with two 
EMI 9124QB bi-alkali photomultiplier tubes to make si- 
multaneous measurements in the Stromgren b and y pass- 
bands. The telescope was programmed to make differ- 
ential brightness measurements of the program star P 



(HD 156846, V = 6.50, B-V = 0.58, Gl V) with respect 
to the two comparison stars CI (HD 157379, V = 6.65, 
B-V = 0.43, F3 IV- V) and C2 (HD 156058, V = 7.70, 
B — V = 0.48, F3 V). To improve the precision of our 
differential magnitudes, we averaged the b and y obser- 
vations to create a (b + y)/2 "passband". The typical 
precision of a single observation in this combined pass- 
band is 0.0015 ma g, as measured from pairs of constant 
stars (|Henrvlll999r i. 

Between 2010, June 01 and July 06, the APT col- 
lected 32 good measurements of the P— CI, P— C2, 
and C2— CI differential magnitudes with standard de- 
viations of 0.00219, 0.00229, and 0.00202 mag, respec- 
tively, slightly higher than the typical 0.0015 mag pre- 
cision. However, HD 156846 and its comparison stars 
are located between —16° and —19° declination, so they 
are observed at higher than average air mass. The ob- 
served scatter of the P-Cl, P-C2, and C2-C1 differ- 
ential magnitudes are all consistent with constant stars. 
Periodogram analyses of the three data sets reveal no 
significant periodicity. We conclude that HD 156846 is 
constant on its rotation timescale. 

6. MONITORING THE TRANSIT WINDOW 

HD 156846 was observed at the Cerro Tololo Inter- 
American Observatory (CTIO) 1.0m telescope using the 
Y4KCam Detector, which is a 4k x 4k CCD with a field of 
view of about 20 arcminutes on the sid43- We monitored 
HD 156846 in the Johnson B band on the nights of 3 Sep 

http:/ /www. astronomy.ohio-state.edu/Y4KCam/ 
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Fig. 3. — Binned photometry of HD 156846 from the night of the 
transit window. The overlaid solid line shows the predicted transit 
signature. 

2009 - during its transit - and 4, 8, 9, 11, and 12 Sep for 
out-of-transit calibration purposes. Due to HD 156846's 
brightness (V = 6.506; B = 7.063), we used a diaphragm 
in the shape of a ring constructed of plywood to decrease 
the effective aperture of the telescope, thereby reducing 
flu x by around 40%, sim ilar to the technique employed 
bv lLopez- Morales! (|2006l ). For the observing strategy and 
ph otometry, our met h ods a re similar to those described 
bv iSouthworth et all ((2009). To ensure that no part of 
the stellar point spread function (PSF) reaches the non- 
linearity levels of the CCD, and to maximize the num- 
ber of counts contained in the PSF, we defocused the 
telescope such that the full-width-half-maximum of the 
PSF subtended 3-8 arcseconds, which produced ~ 10 6 
ADUs (and thus 1 millimag photon noise) per PSF per 
measurement. Observing conditions on the night of the 
transit window (2009 September 3) were not optimal 
and plagued with thin cirrus. We selected the bright- 
est four reference stars in the frame and relative pho- 
tometry was performed using the methods described in 
lEverett fe Howell (l200ll) . 

The new calculated stellar radius presented in this pa- 
per is 2.12i?0, significantly larger than predictions based 
purely upon the spectral type and luminosity class of the 
star. This has three primary effects on the transit pre- 
diction. The first is to increase the transit probability 
to 5.4%. The second is to decrease the predicted tran- 
sit depth to 3 millimags. The third is to increase the 
predicted transit duration to ~ 9 hours. The last two as- 
pects are particularly harmful to attempts at detecting 
a transit since they increase the photometric precision 
requirements and decrease the chances that the window 
can be monitored while the target is observable. 

Figure [3] shows the photometry from the night of the 
transit window. To improve the rms scatter, we binned 
the data into 30 equal time intervals. This improved the 
la scatter from 9 millimags to 7 millimags. Unfortu- 
nately, the poor conditions on that night prevented the 
necessary precision from being achieved since this is still 
a factor of two greater than the predicted transit depth. 
We calculated the pred icted transit signature b ased upon 
the analytic models of IMandel fc Ago! ()2002i ). overplot- 
ted as a solid line in the figure. For a more detailed study 
on the effect of ecc entric orb i ts on transi t lightcurves, w e 
refer the reader to iKippingl (|2008l ) and Kipping] (|2010D . 
Although we see no evidence for a transit in our data, the 
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Fig. 4. — The maximum orbital inclination for a non-transiting 
planet as a function of the argument of periastron for eccentricities 
of 0.0 (solid line), 0.6 (dashed line), and 0.85 (dotted line), plotted 
for periods of both 50.0 days and 359.5 days. The vertical line in- 
dicates the location of the measured periastron argument reported 
here. 

photometric precision is inadequate to rule out such an 
event. For the stellar radius we adopt, the predicted tran- 
sit depth is likely to be q uite robust ag a inst v ariations 



in the planetary radius. iFortnev et al.l (|2007ft showed 



that, for a given planetary composition, planetary radii 
should not vary substantially between orbital radii of 0.1- 
2.0 AU. In order to produce a transit depth comparable 
to the level of precision on the night of the transit win- 
dow, the radius of the companion would need to be > 1.8 
Jupiter radii. Furthermore, the increase in the estimated 
stellar radius caused the predicted end of the transit to 
occur when the lightcurve is heavily influenced by both 
high airmass and cloud contamination. Thus improved 
precision would not have allowed sufficient in and out of 
transit data to be acquired in order to comfortably secure 
the detection. 

If a transit of this planet were to be ruled out, then 
weak constraints on the inclination of the orbit could be 
placed. The magnitude of these constraints as a func- 
tion of periastron argument are shown in Figure |4] which 
shows the maximum orbital inclination for two differ- 
ent periods and three different eccentricities, including 
the period and eccentricity of HD 156846b. A successful 
null-detection would limit the inclination to i < 86.5° 
which would thus place a lower limit on the planetary 
mass of 10.59 Jupiter masses. The period of 50 days is 
shown for comparitive purposes, where one can rapidly 
improve the lower mass limits for the smaller star-planet 
separation. 

7. CONCLUSIONS 

This study has been carried out as part of the 
Transit Ephemeris Refinement and Monitoring Survey 
(TERMS). The purpose of the presented research was to 
improve the orbital parameters of the known exoplanet 
HD 156846b and to monitor the transit window. We 
present new Keck data which, combined with previously 
acquired CORALIE data, refines the orbital parameters 
of the planet. The measurements obtained during pe- 
riastron passage have allowed us to construct an excep- 
tionally accurate transit ephemeris which we present here 
up until the year 2016. The value in a successful tran- 
sit detection would be high for such a long-period planet 
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in an eccentric orbit since it would provide insight into 
the mass-radius relationship for planets in this regime as 
well as allow follow-up characterization to determine the 
radiative time scale and other properties of the atmo- 
sphere. 

The challenge of attempting to monitor the transit win- 
dow is substantial since predicted transit windows are 
so infrequent, the predicted transit duration is relatively 
long, and photometry of bright stars becomes compli- 
cated when comparison stars are few. We present one 
such attempt here where observations were undertaken 
during non-photometric conditions which does not allow 
the transit to be decisively ruled out. A more suitable 
facility to use for such a search is the Microvariability 
and Oscillations of STars (MOST) satellite, such as that 
carried out by (jCroll et al.1 120071 ). From the ground, 
the planned telescopes and instruments of the Las Cum- 
bres Observatory G lobal Telescope (LCOGT) Network 
(jShporer et al.|[20~10h will be ideal for transit monitoring 
due to both their aperture size and longitude coverage. 
The r ecent astrometry work of iReffert &; Quirrenbachl 
appears to indicate that HD 156846b may not be 
in an edge-on orbit, but the results are uncertain enough 
to make this study a worthwhile exercise. We thus en- 
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